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Detection and characterisation of 
coronaviruses in migratory and 
non-migratory Australian wild birds 

Anthony Charmings 1 ' 2 , Tiffanie M. Nelson 1 ' 2 , Jessy Vibin 1 ' 2 , Michelle Wille 3 , Marcel Klaassen 4 & 
Soren Alexandersen 0 1 ' 2 ' 5 

We evaluated the presence of coronaviruses by PCR in 918 Australian wild bird samples collected 
during 2016-17. Coronaviruses were detected in 141 samples (15.3%) from species of ducks, shorebirds 
and herons and from multiple sampling locations. Sequencing of selected positive samples found 
mainly gammacoronaviruses, but also some deltacoronaviruses. The detection rate of coronaviruses 
was improved by using multiple PCR assays, as no single assay could detect all coronavirus positive 
samples. Sequencing of the relatively conserved Orfl PCR amplicons found that Australian duck 
gammacoronaviruses were similarto duck gammacoronaviruses around the world. Some sequenced 
shorebird gammacoronaviruses belonged to Charadriiformes lineages, but others were more closely 
related to duck gammacoronaviruses. Australian duck and heron deltacoronaviruses belonged to 
lineages with other duck and heron deltacoronaviruses, but were almost 20% different in nucleotide 
sequence to other deltacoronavirus sequences available. Deltacoronavirus sequences from shorebirds 
formed a lineage with a deltacoronavirus from a ruddy turnstone detected in the United States. Given 
that Australian duck gammacoronaviruses are highly similarto those found in other regions, and 
Australian ducks rarely come into contact with migratory Palearctic duck species, we hypothesise 
that migratory shorebirds are the important vector for moving wild bird coronaviruses into and out of 
Australia. 


Coronaviruses (CoV) are the causative agents of significant diseases resulting in substantial impact on human 
and animal health. Both severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory 
syndrome coronavirus (MERS-CoV) have caused a significant burden on human health, including a number 
of deaths, and have had significant socioeconomic impacts on the countries in which people were infected ~ 3 . 
Coronavirus infections of livestock, such as porcine epidemic diarrhoea virus (PEDV) and more recently porcine 
deltacoronavirus (PDCoV) in pigs and infectious bronchitis virus (IBV) and turkey coronavirus in poultry, have 
significant impacts on animal health and cause considerable economic costs to producers 4-7 . 

Interspecies spill-over of coronaviruses into new hosts occurs frequently, with SARS-CoV and MERS-CoV 
being the most notable examples of spill-over into humans 8-10 . Bovine coronavirus, canine respiratory corona¬ 
virus, dromedary camel coronavirus and even human coronavirus OC43 all potentially come from the same 
common ancestor, illustrating substantive host flexibility 11-13 . SARS-CoV likely originated in bats while PDCoV 
interestingly is likely to have originated in birds 6,8,9 . Consequently, there is significant interest in assessing wild 
animals for CoV’s. 

Wild birds are ubiquitous and highly mobile potential hosts capable of moving viruses over large distances 
and across geographical and political borders. Wild birds have been implicated in the spread of highly pathogenic 
H5Nx avian influenza viruses 14 and bird migration patterns describe phylogenetic patterns in the matrix gene of 
low pathogenic influenza virus 15 . Coronaviruses have been detected in a range of species of wild birds on all con¬ 
tinents except for Australia and Antarctica 16-21 . It is likely that CoV’s are present in wild birds on all continents. 
However, a recent survey in Australia of 409 birds failed to detect coronaviruses 22 . 
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Studies in which wild bird coronaviruses have been successfully detected, have commonly sampled aquatic 
bird species such as Anseriformes (ducks, geese and swans) and Charadriiformes (gulls, migratory shorebirds 
and waders) 16,18,20,23,24 . This prompted our study to look for coronaviruses in Australian birds that notably belong 
to these two orders. 

Given that the genetic diversity of CoV’s in Australian wild birds was unknown it was decided to use a combi¬ 
nation of PCR assays to maximise the probability of detecting positive individuals. Coronavirus positive samples 
were genetically sequenced and phylogenetic analysis performed to compare these sequences to coronavirus 
sequences from previous surveys of wild and domestic birds across the globe. 

Results 

PCR assay comparison. The S' UTR detected 66/912 or 7.2% (95% Cl: 5.6-8.9%) of samples positive for 
coronavirus, the nested PCR 100/916 samples or 10.9% (95% Cl: 8.9-13.0%) and the modified pancoronavirus 
PCR 109/915 samples or 11.9% (95% Cl: 9.8-14%). By combining the results from the three PCRs, 141/918 or 
15.3% (95% Cl: 13-17.7%) of all samples tested positive for coronaviruses in at least one PCR assay. Nine hundred 
and eight samples were tested with all three PCRs. Ten samples had insufficient RNA to allow testing with all three 
assays. Nine of these samples were negative in the assays with which they were tested. One (Ruddy turnstone 
CoV-9614-2016/11/19-CM/TAS) was positive in the nested PCR, negative in the S' UTR PCR but could not be 
tested with the modified PCR. These ten samples were excluded from the PCR assay comparison. 

When compared to the combined results, the sensitivities for the UTR, nested and modified PCRs were 46.8%, 
70.9% and 77.3% respectively. McNemar’s test revealed no statistical difference in the coronavirus detection abil¬ 
ities of the modified and nested PCRs (\ 2 = 2.38, p = 0.12) but a statistical difference was detected between the 
results of the UTR PCR and the nested PCR — 13.96, p = 0.0002) and the UTR and modified PCR (\ 2 = 21.25, 
p = 4x 10" 6 ). 

From the subset of 66 samples where the Orfl/polymerase PCR amplicon was sequenced and genus of coro¬ 
navirus identified (excluding the one positive pooled sample which contained both a gamma-CoV and delta-CoV 
(see later)), the S' UTR PCR did not detect any delta-CoV’s and only 50% of the gamma-CoV’s (27/54). The 
nested PCR detected 75% (9/12) of the delta-CoV’s and 85% (46/54) of the gamma-CoV’s. The modified PCR 
detected 75% (9/12) of the delta-CoV’s and 87% (47/54) of gamma-CoV’s. No PCR was able to detect all positive 
samples. The S' UTR PCR, despite having the lowest sensitivity, detected 19 samples which were declared negative 
by the other two PCRs. Similarly, both the nested and modified PCRs detected 3 delta-CoV’s the other missed, 
despite targeting exactly the same region of the CoV genome. 

Coronavirus positive samples produced a peak in the melt curve analysis of the modified polymerase real-time 
PCR with delta-CoV positive samples producing a peak between 80.6 °C and 81.8 °C and gamma-CoV positive 
samples producing a peak between 80.4 °C and 85.6 °C. No obvious pattern between melt curve peak tempera¬ 
ture (T m ) or shape was apparent, which could be used to determine the genus of the coronavirus. Samples which 
produced a visible band of the expected size on the agarose gel also produced a strong peak on the melt curve. 
However, some negative samples produced small non-specific peaks in the melt curve within, and slightly above, 
the temperature range in which the positive samples produced peaks. These were small and readily differentiated 
from strong positive samples. However, without running a gel, these peaks would make it difficult to determine 
if a sample was weak positive or a non-specific reactor. When present, these non-specific peaks were associated 
with a -300 bp band on the agarose gel that was readily differentiated from the target band of size 600 bp. During 
optimisation of this PCR on known negative and known positive samples, the non-specific band disappeared with 
forward and reverse primer concentrations of 0.5 pM as opposed to the 1 pM used in this study. This, however, 
had the impact of increasing the C T of positive samples by about 3 (an approximately 8-fold decrease in PCR 
amplification) and it was decided to go for increased PCR sensitivity over the removal of the small non-specific 
product as all positive samples could be confirmed by gel analysis. 

Detection of coronaviruses in Australian wild birds. Coronaviruses were detected in 141 of the 

918 samples tested (15.3%). Positive samples were found in eight out of the 15 bird species tested (Table 1). 
Negative species had very low sample sizes (<30 individuals). Positive species included 4 species of ducks (Order 
Anseriformes): Pacific black duck (18/48), grey teal (19/63), radjah shelduck (1/3) and Australian wood duck 
(1/16); three shorebirds (Order Charadriiformes): curlew sandpiper (8/34), red-necked stint (42/534) and ruddy 
turnstone (47/157); and one species of heron (Order Pelecaniformes): pied heron (5/7). The coronavirus sam¬ 
ple prevalence in Anseriformes was 26.7%, 13.3% in the Charadriiformes and 71.4% in the Pelecaniformes. 
Coronaviruses were only found in apparently healthy birds with none of the dead birds positive for CoV. 
Coronavirus positive birds were found in all states where sampling was conducted. Details of positive samples are 
presented in Supplementary Table SI. 

Of the 141 positive samples, 67 were sequenced to identify the genus of CoV present. Samples for sequencing 
were a subset of samples which were positive in either the nested and/or modified PCRs targeting the polymerase 
gene. Approximately 25-50% of positive samples per species were sequenced, except in the species for which very 
few samples were tested (e.g. the pied herons and the radjah shelduck). In these species with small sample sizes, 
all positive samples were sequenced. Gamma-CoV’s (n = 54) were more frequently identified than delta-CoV’s 
(n = 13). The pooled sample from juvenile Pacific black ducks contained both a gamma- and a delta-CoV. This 
was the only sample in which two coronaviruses were detected. In the order Anseriformes, 22 gamma-CoV’s 
(grey teals (11), Pacific black ducks (10) and radjah shelduck (1)) and one delta-CoV (Pacific black duck) were 
identified. In the order Charadriiformes, 32 gamma-CoV’s (curlew sandpipers (3), red-necked stints (9), ruddy 
turnstones (20)) and 8 delta-CoV’s (curlew sandpiper (1), red-necked stints (2), ruddy turnstones (5)) were 
identified. In order Pelecaniformes, 5 delta-CoV’s were sequenced (pied herons (5)) and no gamma-CoV’s were 
identified. 
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Avian Order/Family 

Species 

Common Name 

coronavirus + ve/Total 

Coronavirus Genus 

Anseriformes, Anatidae 

Anas gracilis 

Grey teal 

19/63 

Gamma (11) 


Anas superciliosa 

Pacific black duck 

18/48 

Gamma (18), Delta 

(D* 


Tadorna radjah 

Radjah shelduck 

1/3 

Gamma (1) 


Chenonetta jubata 

Australian wood duck 

1/16 



Tadorna tadornoides 

Australian shelduck 

0/5 



Anas castanea 

Chestnut teal 

0/11 


Pelecaniformes, Ardeidae 

Ardea picata 

Pied heron 

5/7 

Delta (5) 

Charadriiformes, Scolopacidae 

Caladris ferruginea 

Curlew sandpiper 

8/34 

Gamma (3), Delta (1) 


Caladris ruficolis 

Red-necked stint 

42/534 

Gamma (9), Delta (2) 


Calidris acuminata 

Sharp-tailed sandpiper 

0/3 



Arenaria interpres 

Ruddy turnstone 

47/157 

Gamma (20), Delta (5) 

Procellariiformes 

Ardenna tenuirostris 

Short-tailed shearwater 

0/30 


Acciptriformes 

Milvus migrans 

Black kite 

0/1 


Columbiformes 

Streptopelia chinensis 

Spotted turtle dove 

0/1 


Passeriformes 

Corvus mellori 

Little Raven 

0/4 



Table 1 . The bird species sampled in this study and the proportion of coronavirus positive samples. *One 
pooled sample contained both a gammacoronavirus and deltacoronavirus. 


Phylogenetic analysis of gammacoronaviruses. Gammacoronavirus sequences from this study were 
more related to gamma-CoV’s from ducks and wild birds (Fig. 1) than to IBV-like gamma-CoV’s including IBVs 
from Australia (Supplementary Figure SI). All Australian duck gamma-CoV’s grouped on two branches within 
a predominantly duck lineage of gamma-CoV’s. Bootstrap support for the branches containing the Australian 
sequences was less than 75%. There was no apparent host duck species, temporal or geographical pattern to the 
Australian duck gamma-CoV’s distribution within the phylogenetic tree. 

Gamma-CoV’s from shorebirds grouped on branches with other shorebird gamma-CoV’s sequences obtained 
in this study or from the “Beringia” study of wild birds around the Bering Strait 20 . With the exception of cur¬ 
lew sandpiper gamma-CoV’s, these branches had 75% or greater bootstrap support in the phylogenetic analysis. 
There was sequence diversity in the gamma-CoV sequences in all shorebird species except the curlew sandpipers. 
The majority of the ruddy turnstone gamma-CoV’s grouped with gamma-CoV’s from gulls sampled in 2005 
during the Beringia study One ruddy turnstone gamma-CoV grouped with gamma-CoV’s from rock sandpipers 
from the Beringia study Red-necked stint gamma-CoV’s sequences grouped on their own branch. The three 
curlew sandpiper gamma-CoV sequences were identical, and were most similar to CoV sequences from ducks 
sampled in Sweden during 2011 25 . 

The amino acid ML tree of the gamma-CoV separated the majority of the shorebird gamma-CoV’s and gull 
gamma-CoV’s from the duck-CoV’s (Fig. 2, Supplementary Figure S2). However the three curlew sandpiper 
CoV’s and one ruddy turnstone CoV (Ruddy Turnstone CoV-9564-2016/ll/19-BB/TAS) were identical in amino 
acid sequence to several Australian duck gamma-CoV’s, as well as many domestic and wild duck gamma-CoV 
sequences from Asia, Africa, Europe and North America. 

Phylogenetic analysis of deltacoronaviruses. Deltacoronavirus sequences from this study were spread 
across three lineages in the ML phylogenetic analysis (Fig. 3, Supplementary Figure S3). Delta-CoV’s from 
ruddy turnstones, red-necked stints and curlew sandpipers (all from order Charadriiformes) formed a lineage 
with a delta-CoV sequence from a ruddy turnstone sampled in the United States. The single Australian duck 
delta-CoV belonged to a duck lineage with duck delta-CoV’s detected in Hong Kong but was 18.6% different 
to its nearest relative. The pied heron delta-CoV sequences were all identical, and belonged to a lineage which 
included delta-CoV’s from species from avian orders Pelecaniformes (Herons), Suliformes (Cormorants) and 
Anseriformes (Ducks) from Hong Kong and Cambodia. The pied heron delta-CoV was 17.7% different to the 
next most similar sequence. The genetic diversity in the Charadriiformes delta-CoV lineage was low (1.9-3.2% 
nucleotide differences) compared to the other lineages (duck lineage 18.6-30.9%; heron/cormorant lineage 17.7- 
22.8%). The inter-lineage genetic variation was high between wild bird delta-CoV’s. For example the Australian 
duck delta-CoV was 31.8-34.4% different to the delta-CoV’s within the Heron/Cormorant lineage and 29.6- 
31.8% different to sequences within the Chardriiforme lineage. 

Discussion 

This study is the first to demonstrate the presence of coronaviruses in wild birds in Australia. Coronaviruses 
were distributed across multiple species and geographical locations, with an average coronavirus positive prev¬ 
alence of 15.3%. In bird orders with over 50 individuals sampled, coronaviruses were most frequently detected 
in Anseriformes (26.7%) followed by Charadriiformes (13.3%). Gammacoronaviruses were found in 100% of 
sequenced samples from ducks, while duck deltacoronaviruses were much less common (in only one sample). 
Shorebirds were positive for gamma-CoV’s most commonly (80% of sequenced samples), but delta-CoV’s were 
also present in these species (20% of sequenced samples). Herons were positive for delta-CoV’s exclusively; 
however, with the caveat that only 7 samples were included in this study. While the patterns of predominantly 
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I Duck/Avian CoV - China, Hong Kong 

- KP033077.1 Duck-dominant coronavirus isolate DdCoV/DK/Anhui/1703/2014 

- KT222539.1 Duck-dominant coronavirus isolate DdCoV/DK/jiangsu/FK277/2013 


KJ741882.1 Avian coronavirus strain OH42 Eurasian Teal (2009) 


I Duck CoV - China, Hong Kong 


I Duck/Avian CoV - Hong Kong, Korea & Sweden 

- KT222532.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/HR269/2013 

| Duck CoV - China 

- JN788871.1 Avian coronavirus isolate K561/Northern Shoveler/091223 (2009) 




- JN788786.1 Avian coronavirus isolate J0559/Eurasian Teal/091127 (2009) 


- KP033078 1 Duck-dominant coronavirus isolate DdCoV/DK/Anhui/1656/2014 


I JN788792.1 Avian coronavirus isolate J0579/Eurasian Teal/091127 (2009) 

-^ Avian CoV - Madagascar 


- Ruddy Turnstone CoV-9564-2016/11/19-BB/TAS 

. GU396687.1 Rock sandpiper coronavirus CIR-65824 (2005) 

J GU396689 1 Rock sandpiper coronavirus CIR-65885 (2005) 

I GU396688 1 Rock sandpiper coronavirus CIR-65828 (2005) 

9, JN788808 1 Avian coronavirus isolate J10127/091217 (2009) 

• JN788869.1 Avian coronavirus isolate K554/Northem Shoveler/091223 (2009) 


92 i KM015313.1 Avian coronavirus isolate Mallard Duck/90A88582/133648/111124 (2011) 
I KM015312.1 Avian coronavirus isolate Mallard Duck/90A88582/133583/111123 (2011) 
KM015310 1 Avian coronavirus isolate Mallard Duck/90A88491/131624/111024 (2011) 

KM015309.1 Avian coronavirus isolate Mallard Duck/90A88390/131192/111020 (2011) 

KM015308 1 Avian coronavirus isolate Mallard Duck/90A88309/130049/110916 (2011) 

KM015307.1 Avian coronavirus isolate Mallard Duck/90A88309/129823/110902 (2011) 

KM015306.1 Avian coronavirus isolate Mallard Duck/90A88310/129776/110829 (2011) 

Curlew Sandpiper CoV-9822-2016/12/30-WS/VIC 
Curlew Sandpiper CoV-9819-2016/12/30-WS/VIC 
Curlew Sandpiper CoV-9776-2016/12/28-WS/VIC 
KT882626.1 Avian coronavirus isolate AvCoV/Mallard/70046/2007 
- KY320404.1 Avian coronavirus isolate 937/Mallard Duck/Sweden/2015 


- KY320401.1 Avian coronavirus isolate 485/Mallard Duck/Sweden/2014 

Duck CoV - Sweden 

Red-necked Stint CoV-9772-2016/12/28-WS/VIC 
Red-necked Stint CoV-9752-2016/12/28-WS/VIC 
97 Red-necked Stint CoV-9738-2016/12/28-WS/VIC 
Red-necked Stint CoV-9729-2016/12/28-WS/VIC 
Red-necked Stint CoV-10402-2017/1 /21-BI/VIC 

-Red-necked Stint CoV-9786-2016/12/29-WS/VIC 

, Red-necked Stint CoV-9068-2016/12/30-WS/VIC 
98l Red-necked Stint CoV-10415-2017/1/21-BI/VIC 


- KT222525.1 Duck-dominant coronavirus isolate DdCoV/DK/sichuan/FP276/2013 


Duck/Avian CoV - Hong Kong 


JN788826.1 Avian coronavirus isolate J1362/Northem Shoveler/091223 (2009) 
JN788830.1 Avian coronavirus isolate J1393/Northern Pintail/091230 (2009) 
JN788833.1 Avian coronavirus isolate J1404/Northem Pintail/091230 (2009) 
JN788862.1 Avian coronavirus isolate J1581/Northem Pintail/100112 (2010) 

GU396668.1 Pintail coronavirus PB/AK-10 (2005) 

GU396671.1 Pintail coronavirus PB/AK-25 (2005) 


-KT222530 1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR266/2013 
- KT222513.1 Duck-dominant coronavirus isolate DdCoV/DK/fujian/67/2013 


KP033043.1 Duck-dominant coronavirus isolate DdCoV/DK/Zhejiang/1411/2014 
- KT882628 1 Avian coronavirus isolate AvCoV/Mallard/70125/2007 


KY320402 1 Avian coronavirus isolate 665/Mallard Duck/Sweden/2014 
KY320400.1 Avian coronavirus isolate 271/Mallard Duck/Sweden/2014 

JN788795.1 Avian coronavirus isolate J0589/Eurasian Teal/091127 (2009) 

JN788784.1 Avian coronavirus isolate J0126/Eurasian Teal/091106 (2009) 

JN788783 1 Avian coronavirus isolate J0124/091106 (2009) 

JN788782.1 Avian coronavirus isolate J0122/Eurasian Teal/091106 (2009) 

JN788781.1 Avian coronavirus isolate J0121/Eurasian Teal/091106 (2009) 

JN788780 1 Avian coronavirus isolate JO 109/Eurasian Teal/091106 (2009) 

KT254275.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F11/2014 

- KP033062.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F11/2014 _ 

Pacific Black Duck CoV-10221-2016/6/13-CR/VIC 
Pacific Black Duck CoV-G0001-2016/12/21-LC/VIC 
Grey Teal CoV-10178-2016/6/13-CR/VIC 


Duck/Avian CoV - China, Hong Kong & Sweden 


- Grey Teal CoV-10181 -2016/6/13-CR/VIC 


Pacific Black Duck CoV-9719-2016/12/22-LC/VIC 
Pacific Black Duck CoV-9711-2016/12/21 -LC/VIC 
-Pacific Black Duck CoV-10180-2016/6/13-CR/VIC 


Grey Teal CoV-10229-2016/10/28-MM/NSW 
Grey Teal CoV-10235-2016/10/28-MM/NSW 

Grey Teal CoV-10250-2016/10/29-MM/NSW 


Grey Teal CoV-10243-2016/10/29-MM/NSW 
Pacific Black Duck CoV-9721-2016/12/22-LC/VIC 
Pacific Black Duck CoV-9718-2016/12/22-LC/VIC 

-Grey Teal CoV-10234-2016/10/28-MM/NSW 

- Pacific Black Duck CoV-8907-2016/6/07-LC/VIC 


- Grey Teal CoV-10002-2016/6/11-PV/VIC 

. JN788831.1 Avian coronavirus isolate J 1395/Northern Pintail/091230 (2009) 
M1JN788835.1 Avian coronavirus isolate J 1407/Northern Pintail/091230 (2009) 
m,GU 396673 1 Pintail coronavirus PB/AK-124 (2005) 

IGU396672.1 Pintail coronavirus PB/AK-37 (2005) 

- KJ741871.1 Avian coronavirus strain OH11 Eurasian Teal (2009) 


KJ741878.1 Avian coronavirus strain OH33 Mallard Duck (2009) 

KM093886.1 Avian coronavirus isolate 7V/Red Fody/2011 

J6| - Pacific Black Duck CoV-9710-2016/12/21-LC/VIC 

i- Radjah Shelduck CoV-9515-2016/4/28-HD/NT 

Grey Teal CoV-10214-2016/6/13-CR/VIC 

. Pacific Black Duck CoV-10196-2016/6/13-CR/VIC 

-^ , Grey Teal CoV-10230-2016/10/28-MM/NSW 

Grey Teal CoV-10228-2016/10/28-MM/NSW 


- KF437666.1 Avian coronavirus isolate 12WB180/Pintail/2012 


- KT882616 1 Avian coronavirus isolate AvCoV/Greater Scaup/67699/2007 


KJ741881 1 Avian coronavirus strain OH38 Mallard Duck (2009) 

- m Duck CoV - China & USA 

JN788794.1 Avian coronavirus isolate J0588/Eurasian Wigeon/091127 (2009) 


76* KT222497.1 Duck-dominant coronavirus isolate DdCoV/DK/anhui/FA271/20l3 
* KP033060.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F28/2014 
KT222536.1 Duck-dominant coronavirus isolate DdCoV/DK/jiangsu/FK278/2013 
KT222540.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR215/2013 

■ KP033059.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F31/2014 

«1- KP033048.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F93/2014 

KT222535.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR217/2013 

- KT222514.1 Duck-dominant coronavirus isolate DdCoV/DK/sichuan/HP229/20l3 


-KT222537.1 Duck-dominant coronavirus isolate DdCoV/DK/anhui/FA246/2013 
KP033070.1 Duck-dominant coronavirus isolate DdCoV/DK/Anhui/F82/2014 

- KT222529.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR228/2013 

KP033038.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/H 1373/2014 


KP033039.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/H 1001-2/2014 

KP033056.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/F34/2014 


«9 KM454473.1 Duck coronavirus isolate DK/GD/27/2014 

KP033040.1 Duck-dominant coronavirus isolate DdCoV/DK/Guangdong/H1001-1/2014 

_I KT222534.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR226/2013 

""571 KT222533.1 Duck-dominant coronavirus isolate DdCoV/DK/hubei/FR225/2013 




- GU396674 1 Glaucus-gull coronavirus PB/AK-173 (2005) 


1 Ruddy Turnstono CoV-9694-2016/11 28-BP/S A 
I Ruddy Turnstone CoV-9696-2016/11/28-BP/SA 
Ruddy Turnstone CoV-9624-2016/11/19-CM/TAS 
Ruddy Turnstone CoV-9612-2016/11 /19-CM/TAS 

- GU396682.1 Glaucous-winged gull CIR-66002 (2005) 



931 Ruddy Turnstone CoV-9592-2016/11/18-BB/TAS 
1 Ruddy Turnstone CoV-9566-2016/11/18-BB/TAS 
Ruddy Turnstone CoV-9662-2016/11/27-NV/SA 
Ruddy Turnstone CoV-9628-2016/11/22-DW/TAS 
Ruddy Turnstone CoV-9617-2016/11/19-CM/TAS 
Ruddy Turnstone CoV-9596-2016/11/19-CM/TAS 

GU396679 1 Black-headed gull coronavirus CIR-66183 (2005) 
GU396680.1 Black-headed gull coronavirus CIR-66185 (2005) 
GU396681.1 Black-headed gull coronavirus CIR-66187 (2005) 
GU396683 1 Black-headed gull coronavirus CIR-66144 (2005) 
GU396686.1 Black-headed gull coronavirus CIR-66152 (2005) 
GU396684.1 Black-headed gull coronavirus CIR-66146 (2005) 


0.01 


Figure 1 . Maximum likelihood phylogenetic tree of the 277 nt fragment of the polymerase gene of 
gammacoronaviruses from wild birds including the sequences obtained in this study Bootstrap confidence of 
each branch was calculated from 1000 replicates. The region where each sequence was obtained is indicated 
with colour. Australia (red), Madagascar (teal), China, Hong Kong and Korea (dark blue), Bering Strait (green), 
United States of America (light blue) and Sweden (pink). The sampling location and state of each Australian 
sample is identified with a two letter code and the state: New South Wales (NSW): MM-Moulamein; Northern 
Territory (NT): HD-Humpty Doo; South Australia (SA): BP-Boatswain Point, NV-Nene Valley; Tasmania 
(TAS): BB-Borges Bay, King Island, CM-Central Manuka, King Island, DW-Dripping Wells, King Island, 
Tasmania; Victoria (VIC): BI-Barrallier Island, CR-Carlisle River, LC-Lake Connewarre, PV-Paynesville, WS- 
Werribee South. Some branches have been collapsed if all sequences came from the same region. The complete 
tree is shown in Supplementary Figure SI. 


gammacoronaviruses in ducks and deltacoronaviruses in herons have been described in previous surveys 16,19,20,23 , 
the detection of both gamma- and delta-CoV’s in the same populations of Charadriiformes at the same time is 
unique to this study. 

DNA sequencing of 67 of the 141 PCR positive samples confirmed that the PCRs were in fact amplifying coro¬ 
navirus cDNA, although it was evident that no single PCR assay could detect all coronavirus positive samples. 
The use of multiple PCR assays to screen wild bird samples increased the overall sensitivity for detecting coro- 
naviruses. Compared to the 5' UTR PCR, nested and modified PCRs alone, the detection rate of coronaviruses 
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Figure 2. Maximum likelihood phylogenetic tree of the 92 amino acid sequence of the gammacoronaviruses 
from wild birds. Bootstrap confidence of each branch was calculated from 500 replicates. Some branches have 
been collapsed for clarity The region where each sequence was obtained is indicated with colour. Australia 
(red), Madagascar (teal), China, Hong Kong and Korea (dark blue), Bering Strait (green), United States of 
America (light blue) and Sweden (pink). The sampling location and state of each Australian sample is identified 
with a two letter code and the state: New South Wales (NSW): MM-Moulamein; Northern Territory (NT): HD- 
Humpty Doo; South Australia (SA): BP-Boatswain Point, NV-Nene Valley; Tasmania (TAS): BB-Borges Bay, 
King Island, CM-Central Manuka, King Island, DW-Dripping Wells, King Island, Tasmania; Victoria (VIC): BI- 
Barrallier Island, CR-Carlisle River, LC-Lake Connewarre, PV-Paynesville, WS-Werribee South. Some branches 
have been collapsed if all sequences came from the same region. The complete tree is shown in Supplementary 
Figure S2. 
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Figure 3. Maximum likelihood phylogenetic tree of the 311 nt fragment of the polymerase gene of 
deltacoronaviruses from wild birds, including the sequences obtained in this study Bootstrap confidence of 
each branch was calculated from 1000 replicates. The region where each sequence was obtained is indicated 
with colour. Australia (red), Madagascar (teal), China, Hong Kong and Korea (dark blue), Bering Strait region 
(green), United States of America (light blue) and Sweden (pink). The sampling location and state of each 
Australian sample is identified with a two letter code and the state: New South Wales (NSW): MM-Moulamein; 
Northern Territory (NT): HD-Humpty Doo; South Australia (SA): BP-Boatswain Point, NV-Nene Valley; 
Tasmania (TAS): BB-Borges Bay, King Island, CM-Central Manuka, King Island, DW-Dripping Wells, King 
Island, Tasmania; Victoria (VIC): BI-Barrallier Island, CR-Carlisle River, LC-Lake Connewarre, PV-Paynesville, 
WS-Werribee South. Some branches have been collapsed if all sequences came from the same region. The 
complete tree is shown in Supplementary Figure S3. 


increased 2.1, 1.4 and 1.3 fold respectively when all PCRs were used together. Had only one assay been used in 
this study, a number of coronaviruses in individual species would not have been detected. For example, only the 
modified PCR detected delta-CoV s in the curlew sandpipers and the red-necked stints, while the nested PCR was 
able to detect more delta-CoV’s in ruddy turnstones and pied herons. None of the PCRs described in this study 
detected the deltacoronavirus in the pooled duck sample. Amplicon sequencing from this sample found that the 
PCRs were preferentially amplifying the gammacoronavirus present. However next generation sequencing data 
from this sample (Vibin et al ., 2018 (submitted)) found deltacoronavirus sequences were more abundant than 
those of the gammacoronavirus. Therefore the current pan-coronavirus PCRs may favour detection of gamma- 
coronaviruses over deltacoronaviruses and this may explain why studies to date have found many more gamma- 
coronaviruses in wild birds. It is possible that there are many more genetically diverse coronaviruses circulating in 
wild birds which are currently outside the detection abilities of the current published pancoronavirus PCRs, and 
it will be interesting to see how our understanding of coronavirus diversity improves as the use of non-targeted 
next generation sequencing techniques becomes more common. 

The small fragment of the coronavirus genome sequenced in this study is one of the most conserved regions 
of the viral polymerase gene 24 . It is the most sequenced genomic region of coronaviruses from wild bird hosts, 
with many studies using PCRs targeting this region to detect avian coronaviruses 16,18,24,26 . The fragment analysed 
in this study represents about 1% of the coronavirus genome, and therefore any interpretation of relationships in 
the phylogenetic tree should be made cautiously. With these caveats in mind, no temporal or spatial patterns in 
coronavirus sequences was observed in either the gamma- or deltacoronavirus phylogenetic analyses. However, 
some clustering by host species/order was observed, and this was strongest in the deltacoronaviruses, similar 
to the situation observed by Chu, et al 16 in Hong Kong and Cambodia who observed similar host species line¬ 
ages. There was a Charadriiformes delta-CoV lineage observed which was strongly supported by the bootstrap 
analysis. All seven Australian shorebird delta-CoV’s from three different host species clustered with the single 
previously characterised shorebird delta-CoV from a ruddy turnstone sampled on the East Coast of the United 
States 21 . The one Australian duck deltacoronavirus belonged to a duck lineage of deltacoronaviruses, but itself 
was 18.6% different to its closest relative sequence from a Eurasian teal in Hong Kong 16 . Similarly, the Australian 
heron deltacoronaviruses were within the lineage containing other herons but were genetically different to heron 
sequences seen elsewhere. 

There were Charadriiformes lineages within the gammacoronavirus phylogenetic tree, and the majority of the 
Australian Charadriiform gamma-CoV s grouped on branches with other gamma-CoV s from host species within 
this order. However some Charadriiform gammacoronaviruses were highly similar in nucleotide sequence to 
gammacoronaviruses from ducks, e.g. the curlew sandpiper gammacoronaviruses contained a single nucleotide 
difference to mallard gammacoronavirus sequences identified in Sweden 25 . Even within the Charadriiformes 
lineages, the degree of genetic difference to duck sequences was much less than that seen in the deltacoronavirus 
phylogeny. The phylogeny based on the amino acid sequence of the gammacoronaviruses separated the majority 
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of the duck and Charadriiform gammacoronaviruses into separate branches, however both the curlew sandpiper 
gammacoronaviruses and a single ruddy turnstone gamma-Co V were identical in amino acid homology to duck 
Australian duck gamma-CoV’s as well as duck gamma-CoV’s identified in other investigations 16,18,20,25,27 . These 
results again suggest that migratory Charadriiformes may share gammacoronaviruses with other avian species 16,18 
and are therefore important vectors and reservoir of coronaviruses. Hughes, et al. 23 and Chu, et al lb identified 
similar gammacoronavirus sequences between Charadriiformes and Anseriformes, while Lima, et al. 18 found 
similar gammacoronaviruses in Charadriiformes and Gruiformes in Madagascar. Interestingly, Hughes, et al 23 
identified their similar sequences by a PCR targeting the 3' UTR, while Lima, et al 18 identified similar sequences 
with the same nested PCR used in this investigation 16 . 

Regular spill-over of coronaviruses from Charadriiformes to Anseriformes could explain why many 
Australian duck gammacoronaviruses were similar to the duck gammacoronaviruses from Asia and more distant 
regions. Australian duck species rarely come into direct contact with Palearctic waterfowl 28 . Some degree of geo¬ 
graphical gammacoronavirus lineage formation might therefore have been expected, if this was the only means 
of transmission of duck gamma-CoV’s, due to the high rate of mutation and recombination observed in other 
avian gamma-CoV’s 29,30 . However, such isolation was not observed, which suggests that movement of coronavi¬ 
ruses between these populations of ducks occurs frequently Given that millions of Charadriiformes migrate to 
Australia to spend their non-breeding season in Australia 31 , Australian ducks would be much more likely to come 
into direct contact with Palearctic Charadriiformes than Anseriformes. It is also likely that there are many oppor¬ 
tunities for individuals within the migrating shorebird populations to acquire and transmit coronaviruses to/from 
ducks as they move north-south, as coronaviruses are prevalent in both Charadriiformes and Anseriformes along 
the East Asian-Australasian Flyway 5,20,32 and these species forage and breed in similar habitats 33 . 

This study further highlights the need for greater genetic characterisation of wild bird coronaviruses from all 
regions of the world to better understand the dispersal and host spill-over dynamics of these important pathogens. 
Further work to obtain more sequence data from the Australian wild bird coronaviruses is currently underway 

Materials and Methods 

Sample collection and Storage. Samples tested in this investigation (918) originated from two sources, 
avian influenza A surveillance by Deakin University (912) and a few dead birds submitted through the Wildlife 
Health Surveillance Victoria program (6) at the University of Melbourne’s Faculty of Veterinary and Agricultural 
Science, Werribee, Victoria. Most samples were taken from birds within avian families Anatidae (Anatinae 
and Tadorninae; Order Anseriformes - 146 samples), Scolopacidae (Order Charadriiformes - 728 samples), 
Procellariidae (Order Procellariiformes - 30 samples) and Ardeidae (Order Pelecaniformes - 7 samples). Birds 
were sampled between April 2016 and January 2017 at multiple locations in south eastern Australia and the 
Northern Territory Birds were caught either by hand, cannon netting or baited walk-in traps 34 . Some duck 
samples were obtained from hunter-killed animals. All samples were taken from single individuals and con¬ 
sisted of a combined cloacal and oropharyngeal swab, except for one sample (Pacific Black Duck CoV-GOOOl- 
2016/12/21-LC/VIC), which was a pool of faeces collected from a group of 6 juvenile Pacific black ducks trapped 
together. No birds sampled as part of the influenza surveillance were showing any clinical signs of disease. 
Capture and sampling of birds were carried out in accordance with all relevant guidelines and regulations and 
under animal ethics permits B37-2013 and B43-2016 issued by Deakin University Animal Ethics Committee. 

Oropharyngeal and cloacal swabs from birds belonging to other families (Accipitridae (Order Accipitriformes) 
Columbidae (Order Columbiformes), Corvidae (Order Passeriformes)) were collected from six dead birds. Post 
mortem revealed the cause of death in the four little ravens ( Corvus mellori) and one black kite (. Milvus migrans) 
as physical trauma. Only the spotted turtle dove ( Streptopelia chinensis) showed clinical signs of disease (diar¬ 
rhoea and depression) and had been euthanised by a veterinarian (no post mortem was performed on this bird). 

Oropharyngeal and cloacal swabs were pooled in viral transport media (brain heart infusion [BHI] broth 
based medium [Oxoid] with 0.3 mg/ml penicillin, 5 mg/ml streptomycin, 0.1 mg/ml gentamicin, and 2.5g/ml 
amphotericin B). Swabs were placed in viral transport media at the collection sites, transported at 4°C to the 
laboratory and stored at —80 °C until the nucleic acid extraction was performed. 

Extraction of RNA. Total viral nucleic acids from captured wild bird samples was extracted using the 
Nucleomag Vet Viral DNA/RNA isolation Kit (Scientifix, Australia) and a Kingfisher™ Flex extraction robot 
(Life Technologies) as per the manufacturer’s instructions, with an elution volume of 50 ul. RNA from the six 
samples collected from the dead birds was extracted using the Viral RNA mini-prep kit (Qiagen, Australia) as 
per the manufacturer’s instructions. Extracted RNA was stored at —80 °C. RNA from infectious bronchitis virus 
(IBV) Nl/62 cultured in specified pathogen free chicken eggs (allantoic fluid, passage 7) was extracted using the 
Viral RNA mini-prep kit and used as a PCR positive control in this investigation. The RNA of each sample was 
then used in three separate PCR assays as described below. 

S' Untranslated region (UTR) coronavirus real-time PCR. A real-time TaqMan PCR, originally devel¬ 
oped to detect IBV in commercial poultry 35 , and later used to screen wild birds for the presence of coronavi¬ 
ruses 17 , was adapted into a real-time SYBR Green PCR assay using the Power SYBR Green RNA-to-C T 1 step kit 
(Applied Biosystems, California, USA). The PCR amplified a 142 bp product from the 5' UTR of gammacoro¬ 
naviruses using primers IBV5 / GU391 (5'-GCT TTT GAG CCT AGC GTT-3') and IBV5 / GL533 (5'-GCC ATG 
TTG TCA CTG TCT ATT G-3'). 2 pi of sample RNA was added to 1 x Power SYBR Green RT-PCR Mix, 1 uM 
of each primer, 1 x RT-Enzyme Mix and RNase free water to make a total reaction volume of 10 pi. The reaction 
was conducted in a QuantStudio™ Flex 6 real-time thermal cycler at 48 °C for 30 min, 95 °C for 10 min, 40 cycles 
of 95 °C for 15 sec, 60 °C for 10 sec, 72 °C for 20 sec and a final 72 °C step for 3 min. A melt curve analysis was 
performed immediately post PCR with the reaction conditions of 95 °C for 15 sec, then 60 °C for 1 min followed 
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by a continuous temperature ramp between 60 °C and 95 °C increasing at 0.05 °C/sec. The IBV positive control 
produced a single peak melt curve with a peak melt temperature (T M ) of 80.25 d= 0.3 °C. If a sample produced a 
melt curve with a T M within 0.3 °C of the IBV positive control, 5 pi was then run on a 2% agarose gel. A sample was 
declared positive if it produced a single peak melt curve consistent with the IBV positive control and produced a 
140 bp product on the agarose gel. 

Pancoronavirus polymerase nested PCR. A nested pancoronavirus PCR originally developed to 
detect coronaviruses in bats and later modified and used to detect gammacoronaviruses and deltacoronaviruses 
in wild birds 16 was used to detect coronaviruses in this investigation. The primers used in the initial PCR were: 
Chu-RdRp-Nl-F (5'-GGK TGG GAY TAY CCK AAR TG-3') and Chu-RdRp-Nl-R (5'-TGY TGT SWR CAR 
AAY TCR TG-3 r ) and amplified a 602 bp product from the polymerase gene of coronaviruses. cDNA synthesis 
was performed using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher, Waltham, MA, USA) 
using 1 x random hexamers as per the manufacturer’s protocol. The initial PCR reaction was made with 2 pi of 
cDNA, 1 x AmpliTaq Gold 360 PCR Mix (Applied Biosystems, California USA), 0.5 pM of each primer, 3 mM 
MgCl 2 and nuclease free water in a total volume of 10 pi. The nested PCR used primers Chu-RdRp-N2-F (5'-GGT 
TGG GAC TAT CCT AAG TGT GA-3') and Chu-RdRp-N2-R (5'-CCA TCA TCA GAT AGA ATC ATC AT-3') 
and amplified a 440 bp product within the first PCR. The PCR reaction mixture was made with 1 pi of the first 
PCR reaction, 1 x AmpliTaq Gold 360 PCR Mix, 0.5 pM of each primer, 3 mM MgCl 2 and RNase free water in 
a total volume of 10 pi. PCR reactions were conducted in a thermal cycler at 95 °C for 10 min, then 40 cycles of 
95 °C for 30 sec, 48 °C for 30 sec, 72 °C for 50 sec followed by a final step of 72 °C for 3 min. 5 pi of each PCR was 
then run on a 1.6% agarose gel. A sample was declared positive if a band of the expected size was seen on the gel 
from either PCR. 

Modified pancoronavirus polymerase real-time PCR. A non nested touchdown real-time PCR 
was developed and evaluated against the nested pancoronavirus PCR to see if it could be used as a replace¬ 
ment for the nested assay in future studies. Primers were designed incorporating all current coronavirus poly¬ 
merase gene sequences and looking for the most conserved regions. Primer candidates were then evaluated in 
PerlPrimer 37 . The selected primers, AC-CoV-F (5'-GGT TGG GAT TAT CCW AAR TGT G-3') and AC-CoV-R 
(5'-TGY TGT GAR CAA AAY TCR TG-3'), were less degenerate than Chu-RdRp-Nl-F and Chu-RdRp-Nl-R 
described above and targeted the same conserved sites in the polymerase gene of coronaviruses producing a 
PCR product of 602 bp. cDNA synthesis was performed using the High-Capacity cDNA Reverse Transcription 
Kit (Thermo Fisher, Waltham, MA, USA) as per the manufacturers instructions but using 100 pM of the reverse 
primer AC-CoV-R in place of the 1 x random hexamers. The PCR reaction was performed using 2 pi of cDNA, 
1 x AmpliTaq Gold 360 PCR Mix, 1 pM of each primer, 3 mM MgCl 2 ,2 pM Syto 9 (Thermo Fisher, Waltham, MA, 
USA) and molecular grade water to a final volume of 10 pi. The reaction was conducted in a QuantStudio™ Flex 
6 real-time thermal cycler at 95 °C for 10 min, 45 cycles of 95 °C for 30 sec, annealing temperature for 30 sec, 72 °C 
for 45 sec and a final 72 °C step for 3 min. The touchdown annealing temperature started at 60 °C for 3 cycles, then 
decreased by 2 °C every 3 cycles until 48 °C which was used for the final 30 reaction cycles. 

A melt curve analysis was performed immediately post PCR with reactions subjected to 95 °C for 15 sec, then 
60 °C for 1 min followed by a continuous temperature ramp between 60 °C and 95 °C increasing at 0.05 °C/sec. 
5 pi of the PCR product was then run on a 1.6% agarose gel and a sample was declared positive if a band of the 
expected size was seen. 

PCR product sequencing and phylogenetic analysis. A selection of 67 PCR products from corona- 
virus Orfl/polymerase PCR positive samples were gel purified using the 2% Size Select E-Gel System (Thermo 
Fisher, Waltham, MA, USA) and sequenced using the Big Dye Terminator Cycle v3.1 on a Hitachi 3500x1 genetic 
analyser (Applied Biosystems, Foster City, California, USA). As a result of PCR product sequencing conducted 
during the preliminary PCR testing and optimisation, the pooled duck sample (Pacific Black Duck CoV-GOOOl- 
2016/12/21-LC/VIC) was found to contain both a gamma-CoV and a delta-CoV. The gamma-CoV sequence used 
in this analysis was obtained as described above, however the delta-Co V sequence was obtained from the analysis 
performed by Vibin, et al 38 (submitted). All sequences generated in this study have been deposited in Genbank 
under accession numbers MG764091 to MG764157 and MH090080 and MH090081. 

DNA sequences were aligned in MEGA 7 39 . Each sequence was then trimmed to where both forward and 
reverse reads 100% agreed on the consensus sequence. Consensus sequences between 315 and 515 nt were gener¬ 
ated and queried against the Genbank database using a BLASTn query to identify the genus of coronavirus. Ten 
amplicons which were too short or did not have 100% agreement between the forward and reverse strands were 
excluded from subsequent multiple sequence alignments and phylogenetic analysis. 

A final dataset (57) of 44 gamma-CoV and 13 delta-CoV sequences were aligned with gamma- and delta-CoV 
sequences available on Genbank using the ClustalW algorithm in MEGA 7. The aligned region belonged to the 
Orfl/Polymerase gene between nucleotides 14203 to 14761 of the reference avian infectious bronchitis virus 
sequence (Genbank accession: NC_001451.1). Absence of evidence of recombination was confirmed using 
the Genetic Algorithm Recombination Detection (GARD) method 40 of the HyPhy package hosted on the 
DataMonkey Webserver 41 . Codon specific selection pressure was analysed using the single likelihood ances¬ 
tor counting (SLAC), fixed effects likelihood (FEL), internal fixed effects likelihood (IFEL) and random effects 
likelihood (REL) methods available on the DataMonkey server as previously described 42 . All analyses were 
performed on the optimum model as determined by the model selection algorithm of the server, and using 
Neighbour-Joining phylogenetic trees. A p-value of 0.05 or Bayes factor of 100 were selected for the selection 
pressure analyses, which indicated evidence for negative selection and no evidence of positive selection (data not 
shown). 
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To understand the relationship between the coronavirus sequences identified in this study and those of other 
studies, a maximum likelihood (ML) phylogenetic analysis of the nucleotide sequences with 1000 bootstrap rep¬ 
licates was performed in MEGA 7. The gammacoronavirus phylogenetic tree was created from 277 nt of sequence 
using the Tamura-Nei model with 5 discrete gamma distribution categories 43 . The deltacoronavirus phylogenetic 
tree was created from 311 nt of sequence using the Tamura 3 parameter model 44 with 5 discrete gamma distri¬ 
bution categories with invariant sites. An amino acid maximum likelihood phylogenetic analysis was performed 
on the gammacoronavirus sequence (92 amino acids) using the Jones-Taylor-Thornton model with 5 discrete 
gamma distribution categories and 500 bootstrap replicates. All models were determined to be the best fitting 
models to the sequence data by MEGA 7. 

Statistical and data analysis. A sample was considered positive for coronavirus if it was positive in at least 
one of the PCR assays. The agreement between the pancoronavirus PCRs was evaluated using the McNemar s test 
with Bonferroni’s correction and a significance value (a) of 0.05 5 . 

Data availability. All sequences generated have been deposited in GenBank under accession numbers 
MG764091 - MG764157 and MH090080 and MH090081. Other datasets generated or analysed during the cur¬ 
rent study are available from the corresponding author on reasonable request. 

References 

1. World Bank. People, pathogens and our planet: the economics of one health. (World Bank., Washington, DC 2012). 

2. Lee, S. Costly Lessons From the Middle East Respiratory Syndrome Coronavirus Outbreak in Korea. Journal of Preventive Medicine 
and Public Health 48, 274-276, https://doi.org/10.3961/jpmph.15.064 (2015). 

3. Lee, J.-W. & McKibbin, W. J. In Learning from SARS: Preparing for the Next Disease Outbreak: Workshop Summary The National 
Academies Collection: Reports funded by National Institutes of Health (eds Knobler, S. et al .) (2004). 

4. Cavanagh, D. Coronaviruses in poultry and other birds. Avian Pathol. 34, 439-448, https://doi.org/10.1080/03079450500367682 
(2005). 

5. Ignjatovic, J., Ashton, D. F., Reece, R., Scott, P. & Hooper, P. Pathogenicity of Australian Strains of Avian Infectious Bronchitis Virus. 
/. Comp. Pathol. 126, 115-123, https://doi.org/10.1053/jcpa.2001.0528 (2002). 

6. Ma, Y. et al. Origin, evolution, and virulence of porcine deltacoronaviruses in the United States. American Society for Microbiology 
6, e00064, https://doi.org/10.1128/mBio.00064-15 (2015). 

7. Lee, C. Porcine epidemic diarrhea virus: An emerging and re-emerging epizootic swine virus. Virol J. 12, 193, https://doi. 
org/10.1186/sl2985-015-0421-2 (2015). 

8. Guan, Y. et al. Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China. Science 
302, 276-278, https://doi.org/10.1126/science.1087139 (2003). 

9. Lau, S. K. et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proceedings of the National 
Academy of Science 102, 14040-14045, https://doi.org/10.1073/pnas.0506735102 (2005). 

10. Ferguson, N. M. & Van Kerkhove, M. D. Identification of MERS-CoV in dromedary camels. Lancet Infect. Dis. 14, 93-94, https:// 
doi.org/10.1016/S1473-3099(13)70691-l (2014). 

11. Kin, N. et al. Comparative molecular epidemiology of two closely related coronaviruses, bovine coronavirus (BCoV) and human 
coronavirus OC43 (HCoV-OC43), reveals a different evolutionary pattern. Infect. Genet. Evol. 40, 186-191, https://doi.org/10.1016/j- 
meegid.2016.03.006 (2016). 

12. Vijgen, L. et al. Complete Genomic Sequence of Human Coronavirus OC43: Molecular Clock Analysis Suggests a Relatively Recent 
Zoonotic Coronavirus Transmission Event./. Virol. 79, 1595-1604, https://doi.Org/10.1128/jvi.79.3.1595-1604.2005 (2005). 

13. Lu, S. et al. Discovery of a novel canine respiratory coronavirus support genetic recombination among betacoronavirusl. Virus Res. 
237, 7-13, https://doi.Org/10.1016/j.virusres.2017.05.006 (2017). 

14. Verhagen, J. H., Herfst, S. & Fouchier, R. A. Infectious disease. How a virus travels the world. Science 347, 616-617, https://doi. 
org/10.1126/science.aaa6724 (2015). 

15. Olsen, B. et al. Global patterns of influenza a virus in wild birds. Science 312, 384-388, https://doi.org/10.1126/science.1122438 
(2006). 

16. Chu, D. K. et al. Avian coronavirus in wild aquatic birds. /. Virol.85, 12815-12820, https://doi.org/10.1128/jvi.05838-ll (2011). 

17. Duraes Carvalho, R. et al. Phylogenetic and phylogeographic mapping of the avian coronavirus spike protein-encoding gene in wild 
and synanthropic birds. Virus Res. 201, 101-112, https://doi.Org/10.1016/j.virusres.2015.03.002 (2015). 

18. Lima, F. E. D. et al. Diverse gammacoronaviruses detected in wild birds from Madagascar. European Journal of Wildlife Research 61, 
635-639, https://doi.org/10.1007/sl0344-015-0931-7 (2015). 

19. Wille, M., Muradrasoli, S., Nilsson, A. & Jarhult, J. D. High Prevalence and Putative Lineage Maintenance of Avian Coronaviruses in 
Scandinavian Waterfowl. PLoS One 11, e0150198, https://doi.org/10.1371/journal.pone.0150198 (2016). 

20. Muradrasoli, S. et al. Prevalence and Phylogeny of Coronaviruses in Wild Birds from the Bering Strait Area (Beringia). PLoS One 5, 
el3640, https://doi.org/10.1371/journal.pone.0013640 (2010). 

21. Honkavuori, K. S. et al. Novel Coronavirus and Astrovirus in Delaware Bay Shorebirds. PLoS One 9, e93395, https://doi.org/10.1371/ 
journal.pone.0093395 (2014). 

22. Amery-Gale, J. Investigations into pathogens of potential public health and biosecurity significance in Australian native birds Masters 
of Veterinary Science thesis, The Univeristy of Melbourne, (2016). 

23. Hughes, L. A. et al. Genetically Diverse Coronaviruses in Wild Bird Populations of Northern England. Emerg. Infect. Dis. 15, 
1091-1094, https://doi.org/10.3201/eidl507.090067 (2009). 

24. Muradrasoli, S. et al. Broadly targeted multiprobe QPCR for detection of coronaviruses: Coronavirus is common among mallard 
ducks (Anas platyrhynchos). /. Virol. Methods 159, 277-287, https://doi.Org/10.1016/j.jviromet.2009.04.022 (2009). 

25. Wille, M. et al. Temporal dynamics, diversity, and interplay in three components of the virodiversity of a Mallard population: 
Influenza A virus, avian paramyxovirus and avian coronavirus. Infect. Genet. Evol. 29, 129-137, https://doi.Org/10.1016/j. 
meegid.2014.11.014 (2015). 

26. Duraes-Carvalho, R. et al. Coronaviruses Detected in Brazilian Wild Birds Reveal Close Evolutionary Relationships with Beta- and 
Deltacoronaviruses Isolated From Mammals. /. Mol. Evol., 21, https://doi.org/10.1007/s00239-015-9693-9 (2015). 

27. Zhuang, Q. Y. et al. Genomic Analysis and Surveillance of the Coronavirus Dominant in Ducks in China. PLoS One 10, eO 129256, 
https://doi.org/10.1371/journal.pone.0129256 (2015). 

28. Roshier, D. A., Heinsohn, R., Adcock, G. J., Beerli, P. & Joseph, L. Biogeographic models of gene flow in two waterfowl of the 
Australo-Papuan tropics. Ecol. Evol. 2, 2803-2814, https://doi.org/10.1002/ece3.393 (2012). 

29. Ignjatovic, J., Gould, G. & Sapats, S. Isolation of a variant infectious bronchitis virus in Australia that further illustrates diversity 
among emerging strains. Arch. Virol. 151, 1567-1585, https://doi.org/10.1007/s00705-006-0726-y (2006). 


SCIENTIFIC REPORTS I ( 2018 ) 8:5980 I D0l:10.1038/s41598-018-24407-x 


9 




www.nature.com/scientificreports/ 


30. Quinteros, J. A. et al. Full genome analysis of Australian infectious bronchitis viruses suggests frequent recombination events 
between vaccine strains and multiple phylogenetically distant avian coronaviruses of unknown origin. Vet. Microbiol. 197, 27-38, 
https://doi.Org/10.1016/j.vetmic.2016.ll.003 (2016). 

31. Bamford, M., Watkins, D., Bancroft, W., Tischler, G. & Wahl, J. Migratory shorebirds of the East Asian-Australasian Flyway; 
Population estimates and important sites. (Camberra, Australia, 2008). 

32. Kim, H. R. & Oem, J. K. Surveillance of avian coronaviruses in wild bird populations of Korea. /. Wildl. Dis. 50, 964-968, https://doi. 
org/10.7589/2013-11-298 (2014). 

33. Pizzey, G. & Knight, F. A Field guid to the Birds of Australia. 8th Edition edn, (Harper Collins Publishers, 2007). 

34. Haynes, L. et al. Australian surveillance for avian influenza viruses in wild birds between July 2005 and June 2007. Aust. Vet. J. 87, 
266-272, https://doi.Org/10.llll/j.1751-0813.2009.00446.x (2009). 

35. Callison, S. A. et al. Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis 
virus from infected chickens. /. Virol. Methods 138, 60-65, https://doi.Org/10.1016/j.jviromet.2006.07.018 (2006). 

36. Poon, L. L. M. et al. Identification of a Novel Coronavirus in Bats. /. Virol. 79, 2001-2009, https://doi.org/10.1128/jvi.79.4.2001- 
2009.2005 (2005). 

37. Marshall, O. J. PerlPrimer: cross-platform, graphical primer design for standard, bisulphite and real-time PCR. Bioinformatics 20, 
2471 -2472, https://doi.org/10.1093/bioinformatics/bth254 (2004). 

38. Vibin, J. et al. Detection and characterisation of viruses in Australian wild bird’s faeces: Metagenomics method optimisation and 
initial results. (Submitted) (2018). 

39. Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. 
Evol. 33, 1870-1874, https://doi.org/10.1093/molbev/msw054 (2016). 

40. Kosakovsky Pond, S. L., Posada, D., Gravenor, M. B., Woelk, C. H. & Frost, S. D. W. GARD: a genetic algorithm for recombination 
detection. Bioinformatics 22, 3096-3098, https://doi.org/10.1093/bioinformatics/btl474 (2006). 

41. Delport, W., Poon, A. F., Frost, S. D. & Kosakovsky Pond, S. L. Datamonkey 2010: a suite of phylogenetic analysis tools for 
evolutionary biology. Bioinformatics 26, 2455-2457, https://doi.org/10.1093/bioinformatics/btq429 (2010). 

42. Alexandersen, S., Nelson, T. M., Hodge, J. & Druce, J. Evolutionary and network analysis of virus sequences from infants infected 
with an Australian recombinant strain of human parechovirus type 3. Sci. Rep. 7, 3861, https://doi.org/10.1038/s41598-017-04145-2 
(2017). 

43. Tamura, K. & Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans 
and chimpanzees. Mol. Biol. Evol. 10, 512-526 (1993). 

44. Tamura, K. Estimation of the Number of Nucleotide Substitutions When There Are Strong Transition-Transversion and 
G + C-Content Biases. Mol. Biol. Evol. 9, 678-687 (1992). 

45. Petrie, A. & Watson, P. Statistics for Veterinary and Animal Science. 3rd Edition edn, (John Wiley & Sons, 2013). 

Acknowledgements 

The work was supported by funding from Deakin University, Barwon EEealth and CSIRO and from NEEMRC 
Equipment Grant GNT9000413 to S.A. The Melbourne WEEO Collaborating Centre for Reference and Research 
on Influenza is funded by the Australian Department of Eiealth. The sampling by Deakin University importantly 
relied on the assistance of a large number of volunteers within the Victorian Wader Study Group (sampling 
shorebirds) and Field & Game Australia (ducks). Eiumpty Doo Barramundi (NT), Phillip Island Nature Parks 
(VIC) and Trent Leen (VIC) are thanked for providing access to their land and facilities and providing logistic 
support during field sampling. The sampling was financially supported by ARC discovery grants DP 130101935 
and DP160102146 to M.K. Samples from dead birds were kindly provided by W. Steele from Ecology Balance, 
Victoria, Australia and P. Whiteley from the Faculty of Veterinary and Agricultural Sciences, The University of 
Melbourne, Victoria, Australia. Infectious Bronchitis Virus (IBV) Nl/62 for the PCR positive control was kindly 
provided by D. Williams, S. Sapats and B. Vanderheide at the Australian Animal Eiealth Laboratory, CSIRO, 
Geelong, Victoria, Australia. 

Author Contributions 

A.C. designed the modified PCR and performed the PCRs and DNA sequencing. M.W. collected samples and 
extracted viral RNA. J.V. performed next generation sequencing on the pooled sample and provided the duck 
deltacoronavirus sequence. T.N. performed initial testing of pancoronavirus PCRs. M.K. collected samples and 
coordinated the bird sampling. S.A. coordinated all work done at GCEID. A.C. wrote initial manuscript with 
input from S.A. and later versions were based on input and suggestions from all. 

Additional Information 

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-24407-x. 
Competing Interests: The authors declare no competing interests. 

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and 
institutional affiliations. 

Open Access This article is licensed under a Creative Commons Attribution 4.0 International 
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or 
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre¬ 
ative Commons license, and indicate if changes were made. The images or other third party material in this 
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the 
material. If material is not included in the article’s Creative Commons license and your intended use is not per¬ 
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the 
copyright holder. To view a copy of this license, visit http://creativecommons.Org/licenses/by/4.0/. 

© The Author(s) 2018 



SCIENTIFIC REPORTS I ( 2018 ) 8:5980 I D0l:10.1038/s41598-018-24407-x 


10 






